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摘要 : 为 探讨 第 四 纪 冰 期 气候 波动 对 青藏 高 原 多 年 生 草 本 植物 青海 当归 地 理 分 布 格局 形成 的 影响 ， 检 测 了 


青海 当归 16 个 居 群 147 个 个 体 的 核糖 体内 转录 间隔 区 (ITS) 的 序列 变异 ， 共 发 现 了 7 个 单 倍 型 。 结 果 表 























明 除 松潘 居 群 外 其 余 15 个 居 群 都 拥有 单 倍 型 H1， 青 藏 高 原 边缘 的 居 群 拥有 6 种 ITS 单 倍 型 ， 其 中 5 种 为 
私有 单 倍 型 。 谱 系 关系 分 析 显 示 单 倍 型 Hl 为 较 古老 的 单 倍 型 ，AMOVA 分 析 表 明 ， 居 群 间 的 变异 占 了 总 

















变异 量 的 67% ， 居 群 遗传 分 化 显著 ， 











并 且 存 在 较 明显 的 谱系 地 理 结构 (Ns,>Gs，P<0.05)。 我们 认为 青 









































海 当归 在 第 四 纪 冰 期 最 大 一 次 冰期 来 临时 退缩 到 青藏 高 原 边 缘 的 避难 所 ， 冰 期 后 在 扩张 过 程 中 经 历 了 英 基 





者 效应 形成 了 如 今 的 分 布 格局 。 
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Phylogeography of Angelica nitida (Apiaceae ) ， 


an endemic to the Qinghai-Tibet Plateau 
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Abstract; In this study, we examined nuclear ribosomal DNA internal transcribed spacer (ITS) sequence variation 
of Angelica nitida, a perennial plant species that is endemic to alpine scree in the Qinghai-Tibetan Plateau. Seven 
ITS haplotypes were detected among 147 individuals sampled across 16 populations. All populations except one 
(population SP) contained haplotype ( H1). However, analysis of molecular variance (AMOVA) showed that ITS 
variation was mainly present among populations (67% ) and that significant phylogeographical structure was evident 
(Nsr>Gsr, P<0.05). We inferred from our results that the current distribution pattern of A. nitida resulted from range 
retreats and expansions during the Quaternary. During Quaternary glaciations A. nitida retreated to refugia at the edge 
of the Qinghai-Tibet Plateau , but then expanded its range during interglacials and the current post-glacial period. 
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植物 对 气候 和 环境 变化 的 响应 表现 为 改变 其 。” 第 四 纪 冰 期 和 间 冰 期 的 交互 作用 下 ,青藏 高 原 植 
地 理 分 布 区 或 适应 性 进化 (Hallam, 1994; Etter- ”被 的 变化 不 仅 表现 为 植被 类 型 的 演 奉 ， 还 表现 为 
son 和 Shaw, 2001 ) 。 而 前 人 的 研究 则 表明 对 现代 ”物种 的 迁移 和 分 布 区 的 变化 (Hewitt, 2001)， 很 
生物 类 群 的 分 布 格局 影响 最 大 的 是 第 四 纪 ( 约 2 多 类 和 群 经 历 了 迁移 -扩张 -再 迁移 -再 扩张 的 过 程 ， 
百 万 年 前 ) 以 来 的 气候 变化 (Hewitt, 1996), Æ ”而 且 这 个 过 程 极 大 地 影响 了 物种 的 现代 分 布 格局 
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和 遗传 结构 (Abbott 等 , 2000; Petit 和 Grivet, 2002) 。 

被 称 为 “世界 屋 疹 ”和 “世界 第 三 极 ” 的 
青藏 高 原 ， 也 称 为 全 球 气候 变化 的 “驱动 机 ” 
和 “放大 器 ”( 潘 保 田 和 李 吉 钧 ,1996) ， 是 世界 
地 质 史上 面积 最 大 、 海 拔 最 高 的 高 原 ( 施 雅 风 
等 , 1998)。 最 早 且 较为 全 面 系统 的 阐述 该 地 区 
谱系 历史 形成 的 研究 是 利用 叶绿体 DNA 对 高 原 
台面 和 边缘 的 祁 连 圆 柏 的 研究 (Zhang 等 ， 
2005 ) 。 近 年 来 有 关 青 藏 高 原 及 附近 地 区 气候 变 
迁 与 高 山 物 种 的 谱系 地 理 研 究 发 展 十 分 迅速 ， 已 
有 的 研究 发 现 ， 植 物 的 地 理 分 布 格局 往往 反映 出 
他 们 对 过 去 冰期 气候 震荡 的 响应 ( 张 蔚 等 ， 
2005; Meng 等 , 2007; Chen 等 , 2008 ; Yang 等 ， 
2008; Wang 等 , 2009a, b; Wu 等 , 2010)。 青 藏 
高 原 的 隆 升 以 及 第 四 纪 冰 期 不 仅 对 我 国 西部 地 区 
及 东亚 的 气候 、 自 然 环 境 和 生态 造成 了 极 大 的 影 
响 (Zhang, 1983), ， 而 且 对 全 球 气 候 和 环境 变 
化 等 也 产生 了 深刻 作用 (Raymo 和 Ruddiman, 
1992) 。 因 此 研究 该 区 域内 物种 的 遗传 结构 和 种 
群 演化 历史 有 助 于 理解 该 地 区 的 地 质 和 气候 的 演 
变 历 史 (Hewitt, 2000) 。 

近年 来 核 基 因 组 的 遗传 信息 也 不 断 地 用 于 推 
测 物 种 谱系 格局 的 研究 (PSL, 2009; 王 晓 雄 
等 , 2011 ) 。nrDNA 转录 间隔 区 (Internal tran- 
scribed spacer, ITS) 具有 进化 速率 较 快 、 测 序 难 

















度 较 小 的 特点 ， 是 研究 较 低 分 类 阶 元 系统 发 育 关 
系 应 用 最 为 广泛 的 核 基 因 分 子 标记 之 一 ( Bald- 
win 等 , 1995 ) ， 近 年 来 已 应 用 于 推测 物种 谱系 历 
史 的 研究 。 

青藏 高 原 特有 种 青海 当归 (Angelica nitida 
Wolff) 是 典型 的 靠 种 子 传播 扩散 的 多 年 生 草 本 ， 
主要 分 布 于 四 川西 部 、 甘 肃 南 部 以 及 青海 东部 ， 
生长 于 海拔 2 600 ~4 000 m 的 高 山 草 旬 和 灌 从 中 
(Sheh 等 , 2005)。 前 人 人 研究 表明 青海 当归 在 第 三 
纪 上 新 世 才 开始 出 现 ， 其 所 在 的 青藏 世系 的 演化 过 
程 与 青藏 高 原 的 近代 抬升 有 关 (Liao 等 , 2012), 
本 研究 由 在 通过 对 青海 当归 147 个 个 体 的 ITS 基 
因 片 段 的 测序 分 析 来 揭示 该 植物 单 倍 型 的 地 理 分 
布 式样 、 居 和 群 遗 传 结构 等 ， 进 而 推测 其 在 第 四 纪 
冰期 时 可 能 的 避难 所 和 冰期 后 的 迁移 路 线 。 














1 材料 和 方法 
1.1 实验 材料 

研究 材料 采 自 四 川西 部 、 青 海 东 部 和 甘肃 南部 ， 履 
盖 了 该 种 的 所 有 分 布 区 。 按 居 群 采样 标准 对 青海 当归 进 
行 取 样 ， 即 每 个 居 群 随机 采样 4~11 个 个 体 ， 且 每 个 居 
群 内 采样 个 体 之 间 至 少 间隔 20 m， 共 采集 了 青海 当归 16 
个 居 群 147 个 个 体 (R1, 图 1)。 采 集 植物 生长 良好 的 
幼 叶 ， 用 硅胶 快速 干燥 后 带 回 实验 室 。 和 凭证 标本 存 于 四 
川 大 学 标本 馆 (SZ) 。 









































R1 青海 当归 16 个 居 群 材料 来 源 及 各 居 群 单 倍 型 的 组 成 和 个 体 数 


Table ] Origin, sample size, haplotypes and number of individuals per haplotype for the 16 sampled population of Angelica nitida 















































居 群 采样 地 Populations 样本 数 经 度 4 度 海拔 单 倍 型 及 相应 个 体 数 
No. Longitude Latitude Altitude/m Haplotypes (No. of individul) 

URL (Gansu, Mixian) MX 10 104°33' 34°19" 2 600 H1(10) 

碌 曲 (Gansu, Luqu) LQ 7 102°4" 33°59’ 3 875 H1(6) H3(1) 
夏河 (Gansu, Xiahe) XH 10 102°31' 35°24" 3155 H1 (10) 

榆 中 (Gansu, Yuzhong) YZ 11 103°57' 35°46’ 3 125 H1(3) H2(8) 
IT (Qinghai, Menyuan) MY 10 101°23' 37°28" 3 150 H1(10) 

大 通 (Qinghai, Datong) DT 11 101°34' 37°6' 3 100 H1(11) 

泽 库 (Qinghai, Zeku) ZK 10 101°27' 3591" 3512 H1(10) 

河南 (Qinghai, Henan) HN 4 101°36' 34°44’ 3 640 H1(4) 

同 德 (Qinghai, Tongde) TD 10 100°47' 34°43’ 3 425 H1(9) 

班 玛 (Qinghai, Banma) BM 7 100°51' 32°32! 3 800 H1(7) 

JAH (Qinghai, Jiuzhi) JZ 9 100°56’ 33°31’ 3 800 H1(8) H3(1) 
红 原 (Sichuan, Honyuan) HY 9 102°32" 32°47" 3 850 H1(7) H4(2) 
瓦 切 (Sichuan, Waqie) WQ 10 102°37' 33°7' 3 600 H1(8) H3(2) 
唐 克 (Sichuan, Tangke) TK 9 102°28" 33°24’ 3 400 H1(6) H3(2) H5(1) 
红星 (Sichuan, Hongxing) HX 10 102°44" 34°5' 3 200 H1(10) 

松潘 (Sichuan, Songpan) SP 11 104°7' 32°34’ 4 000 H6(10) H7(1) 
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图 1 青海 当归 ITS 单 倍 型 在 各 
Fig. 1 


1.2 DNA 提取 和 PCR 扩 增 

采用 改良 的 CTAB YE (Doyle 和 Doyle, 1987) 提取 
EA DNA, FH ITS4 和 ITS5 通用 引物 (White 等 , 1990) 对 
nrDNA 的 ITS 片段 进行 PCR P34, PCR 反应 体系 在 30 
AL 体系 中 进行 ， 包 括 模 板 DNA10-40 ng, 3 pL PCR 反 
应 缓冲 液 ，3 pL MgCl, (25 mmol .1L-) ，3 pL dNTP 混合 
液 (10 mmol -L') ， 正 反 向 引物 各 1.5 uL (10 uL mol - 
L''), TaqDNA RAM 0.3 单位 (Tiangen, PE), S 
程序 如 下 : 94 CEFE 4 min， 随 后 94 °C 45s, 54 
45s, 72C 1min, JEI 35 次 ， 最 后 在 72 % 条 件 下 延伸 
7 min。 扩 增 得 到 的 PCR 产物 在 1% 的 琼脂 糖 胶 上 电泳 检 
测 ， 送 上 海 英俊 测序 公司 测序 。 
1.3 数据 分 析 

用 CLUSTAL X (Jeanmougin 等 , 1998) 软件 对 测定 
的 147 个 个 体 的 ITS 全 序列 进行 对 位 排列 ， 并 加 以 手工 
校对 。 用 DNASP5.0 (Librado 和 Rozas, 2009) 确定 单 倍 
型 。 根据 每 个 居 群 内 不 同 单 倍 型 的 个 体 数目 和 总 的 个 体 
数 来 计算 单 倍 型 频率 ， 依 照 Pons 和 Petit (1996) 所 描述 
的 方法 ， 用 PERMUT (available at http://www. pierroton. 
inra. fr/genetics/labo/Software/) 软件 计算 居 群 内 的 遗传 
多 样 性 fH 和 居 群 间 平 均 遗 传 多 样 性 HPA J RP EG 
化 系数 Gs 和 Ns。 居 群 间 的 Gs 和 Ns, 值 用 1 000 次 置换 






























































居 群 中 的 分 布 图 ， 居 和 群 简称 见 表 1 


Geographical distributions of haplotypes detected in A. nitida populations. Population names correspond to those listed in Table 1 











进行 检验 。 应 用 ARLEQUIN 软件 包 version 3.1 ( Exeoffier 
等 , 2006) 中 的 AMOVA (Analysis of Molecular Variance) 
计算 居 群 内 和 居 群 间 的 遗传 变异 水 平 ， 以 及 单 倍 型 分 布 
HY Fyfe. (Weir FU Cockerham, 1984) ， 以 进一步 分 析 居 群 
的 分 化 程度 (1 000 次 置换 检验 )。 应 用 Network 4. 6.0.1 
(Polzin 和 Daneshmand, 2003) 对 青海 当归 检测 到 的 16 
个 ITS 单 倍 型 进行 基于 简约 性 原则 的 中 间 连 接 网 法 分 析 
(Median joining networks) (Bandelt 等 , 1999) ， 并 构建 单 
倍 型 网 状 图 。 以 峨眉 当归 (Angelica omeiensis) 和 长 序 当 
日 (Angelica longopes) 为 外 类 群 ， 用 软件 BEAST 1.6.1 
(Drummond 和 Rambaut, 2007) 构建 各 单 倍 型 的 系统 发 
育 关 系 ， 模 型 为 GTR+I+G。 物 种 水 平 上 居 群 间 的 基因 流 
(Nm) 和 Tajima’s D (Tajima, 1989) 及 Fu's F, (Fu, 
1997) 两 种 无 限 突 变 位 点 模型 的 中 性 检验 方法 检测 及 失 
配 分 析 (Rogers 和 Harpending, 1992) 分 别 在 软件 AR- 
LEQUIN 和 DNASP5. 0 中 完成 。 失 配 分 布 分 析 主 要 用 不 
同 的 单 倍 型 之 间 的 差异 分 布 来 评估 区 域 范 围 内 种 群 近期 
是 否 经 历 过 扩张 ， 当 失 配 分 布 呈 现 为 高 低 不 齐 或 几何 分 
布 时 ,说 明 种 群 没有 经 历 扩 张 ， 当 失 配 分 布 呈 现 单 峰 钟 
形 分 布 曲线 ， 说 明 种 群 近期 经 历 了 扩张 。 我 们 还 通过 观 
察 分 布 和 期 望 分 布 之 间 的 偏离 总 方差 (sum of squared de- 
viation, SSD) 来 检验 种 群 扩张 ， 当 观察 偏离 总 方差 与 期 
望 偏离 总 方差 没有 显著 性 差异 时 ， 即 P>0.05, 符合 种 群 
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扩张 模型 假设 (Rogers,，1995) 。 同 时 使 用 Happending’s 
Ragedness index (RAG) 检验 失 配 分 布 曲线 的 平滑 度 。 

如 果 种 群 符合 种 群 扩张 模型 ， 则 利用 公式 : r=2ut 
来 计算 种 群 扩张 的 近似 世代 数 t (Rogers 和 Hapending， 
1992), FER u 是 每 条 序列 每 代 的 突变 率 ， 通过 u=2 pkg 
计算 得 到 u，k 是 每 个 碱 基 的 突变 率 ，k 是 被 分 析 序 列 的 
平均 碱 基数 ，g 代表 物种 每 个 世代 所 需 时 间 。 在 本 研究 
中 ,本文 的 ITS 序列 的 平均 碱 基数 为 600 bp， 每 个 世代 
的 时 间 保 守 估 计 为 3 年 (根据 我 们 在 野外 的 观察 和 栽 
培 ) ， 该 种 ITS 基因 的 平均 进化 速率 利用 Liao 等 (2012) 
的 分 子 钟 数据 用 Beast 和 Tracer 软件 进行 反 推 ,得 到 的 
该 种 的 ITS 基因 的 平均 进化 速率 为 3.5x10”s/s/y (95% 
的 置信 区 间 为 2.02-5.52x10”s/s/y)。 
























































2 结果 

对 青海 当归 的 16 个 居 群 147 个 体 的 nrITS 进 
行 了 PCR 扩 增 片段 的 全 序列 测定 。 对 位 排列 后 ， 
序列 全 长 604 bp, ， 检 测 出 6 处 变异 位 点 ( 表 2)， 
由 软件 DNASP5.0 确定 了 7 个 ITS 单 倍 型 ( Hl- 
H7) ， 单 倍 型 在 各 居 群 的 分 布 见 图 1。 在 得 到 的 
7 个 单 倍 型 中 ， 单 倍 型 Hl 分布 最 广 ， 所 占 的 比 
例 最 大 (80.95% ) ， 除 四 川 松潘 居 群 外 所 有 的 
青海 当归 居 群 都 有 单 倍 型 H1， 且 在 每 个 居 群 中 
都 为 主要 单 倍 型 ; 单 倍 型 H3 分 布 于 碌 曲 ， 久 
治 ， 唐 克 ， 瓦 切 这 4 个 居 群 中 ; 其 他 5 种 单 倍 型 
都 只 分 布 在 单个 种 群 中 ， 单 倍 型 HH，HS5 ，H6 ， 
H7 则 为 青藏 高 原 东南 部 居 群 的 私有 单 倍 型 ( 表 
1， 图 1) 。 根 据 居 群 在 青藏 高 原 范 围 的 具体 地 理 
位 置 ， 我们 将 青海 当归 的 16 个 居 群 分 为 青藏 高 
原 边缘 地 区 ORE, PR, A, MI, R, 



































德 ) 不 同 的 地 理 居 群 ， 青 藏 高 原 台 面 不 具有 任 
何 的 私有 单 倍 型 。 在 所 研究 的 16 个 居 群 中 ， 青 
藏 高 原 东南 部 横断 山区 的 居 群 ( 红 原 ， 瓦 切 ， 
唐 克 ， 红 星 ， 松 潘 ) 拥有 6 种 单 倍 型 ， 其 中 唐 
克 居 群 在 所 有 居 群 中 的 遗传 多 样 度 最 高 ( 表 1 ， 
图 1)。 

以 峨眉 当归 (Angelica omeiensis ) 和 长 序 当 
归 (Angelica longipes ) VE 为 外 类 和 群 ， 利用 BEAST 
的 Bayes 方法 构建 的 青海 当归 单 倍 型 的 系统 树 如 
图 2A， 结 果 表 明 青 海 当 归 的 7 个 单 倍 型 形成 3 
个 分 支 ， 青 藏 高 原 东 南部 的 红 原 和 唐 克 居 群 的 两 
个 私有 单 倍 型 H4 和 HS 形成 一 个 分 支 ， 位 于 系 
统 树 的 基部 ; 松潘 居 群 的 2 个 私有 单 倍 型 H6 和 
H7 分 支 与 其 他 剩余 单 倍 型 (H1, H2 和 了 H3) 分 
支 构成 姐妹 群 。 
AMOVA 分 析 结 果 ( 表 3) 表明 ， 青 海 当 归 
在 居 群 间 的 变异 (67%) 高 于 居 群 内 的 变异 
(33% ) ， 居 和 群 间 存 在 显著 遗传 分 化 (Fo, =0. 67, 
P<0.001) ， 居 群 间 的 平均 基因 流 Nu 为 0. 245。 
通过 PERMUT2 程序 计算 得 出 ， 居 群 内 平均 遗传 
多 样 性 五 (+se)、 总 的 遗传 多 样 性 H, (se), 
居 群 间 遗 传 分 化 G (+tse) 和 Ng (+tse) 分 别 
H: 0.152(+0.0489), 0.317( +0. 1146), 0.521 
(+0.0287)，0.638 (+0.1681),U itt m 
AR, North BRF Coy (P<0.05) ， 表 明 青 海 当 归 
亲缘 关系 相近 的 单 倍 型 在 相同 或 相近 的 居 群 中 发 
生 ， 表 现 出 明显 的 谱系 地 理 结构 。 通 过 Network 
软件 构建 的 网 状 进化 树 显示 7 种 单 倍 型 之 间 的 系 
统 关 系 见 图 2B。 单 倍 型 Hl 处 于 网 络 树 的 中 心 ， 

















































































































红星 ， 榆 中 ,大 通 ， 门 源 ) 和 青藏 高 原 人 台面 
CZE, WW, SH, A, W, RI, E 


可 能 是 较为 古老 的 单 倍 型 。 
ITS 单 倍 型 碱 基 变 异 的 失 配 分 布 分 析 表 明 ， 


表 2 青海 当归 7 种 单 倍 型 的 变异 位 点 及 个 体 数 


Table 2 Variable sites of the aligned sequences of ITS fragment in the 7 haplotypes of A. nitida 











单 倍 型 Haplotypes 变异 位 点 Variable sites 个 i ik 数 
55 255 257 402 551 587 No. of individual 
H1 A G G c c T iis 
H2 A G T G c T 8 
H3 i T G G c T 6 
H4 G G G G c T 2 
H5 G G G A c T i 
H6 A G G c A č ið 
H7 A G G G c c i 
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单 倍 型 之 间 成 对 碱 基 差 异 出 现 的 观察 频率 与 期 望 不 显著 (P>0.05) (484), 但 AMOVA 分 析 结 

频率 趋势 相似 ， 呈 明显 的 单 峰 分 布 (图 3), 失 ”表明 有 效 种 群 大 小 变化 范围 (pm to % ) 非常 高 
配 分 布 的 偏离 总 方差 (SSD) 和 平滑 度 (RAG) ( 表 4)， 所 以 本 文 认为 青海 当归 种 群 经 历 了 扩 
指数 (P>0.05) 都 支持 青海 当归 居 群 经 历 过 快速 ” 张 。 利 用 当归 属 ITS 的 进化 速率 推算 得 到 的 青海 
的 居 群 扩张 ， 中 性 检验 (Tajima’s D、Fu’s Fo) 当归 居 群 扩张 的 时 间 为 0. 138 百 万 年 前 。 
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A. longipes 





A. omeiensis 





图 2 A. 基于 青海 当归 7 个 ITS 单 倍 型 构建 的 系统 发 育 关系 树 ， 分 支 上 的 数字 分 别 表示 贝 叶 斯 后 验 概率 (BPP) ; 
B. 青海 当归 7 个 ITS 单 倍 型 的 网 状 关系 树 ， 树 中 圆圈 的 大 小 代表 单 倍 型 频率 


Fig.2 A. Phylogenetic relationships of haplotypes in A. nitida; Values below branches (% ) represent the Bayesian PP values; 














B. Network of haplotypes; sizes of circles in the network are proportional to the observed frequencies of haplotypes 


R3 青海 当归 16 个 居 群 的 分 子 遗 传 变异 分 析 
Table 3 Hierarchical analyses of molecular variance (AMOVA) based on ITS of A. nitida 
































RE 变 | Prey 
变异 来 源 总 方差 变异 成 分 in 国定 指数 
M PAE, Degree of a i Percentage of 
Source of variation Sum of squares Variance components 的 全 Fixation index (Fsr) 
freedom variation/ % 
居 群 间 Among populations 15 25.422 0.1755 67.00 
EREN Within populations 132 11.326 0.0865 33.00 Fsr: 0.67 
总 计 Total 147 36.748 0.2620 
Exp 
----@-- Obs 














Pairwise Differences 


图 3 青海 当归 所 有 个 体 nrDNA ITS 片段 数据 的 失 配 分 析 
Fig.3 Mismatch distributions of ITS sequences of A. nitida. The dashed line shows observed values, whereas the solid line 


represents expected values under a model of sudden (stepwise) population expansion 
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表 4 青海 当归 种 群 扩张 的 中 性 检测 
Table 4 Estimates of neutrality tests and expansion time of A. nitida 
SSD RAG Tajima’s D Fu’s Fs Expansion 
Group T % 0, i 
(P value) (P value) (P value) (P value) time( mya) 
Ton 3 5 0 99999 0.00431 0. 24044 -1.12977 -3. 35349 0.138 
a ; (0.5412) (0. 5268 ) (0. 1288) (0. 06160) i 


T: 种 群 群 爆发 以 后 经 过 几 代 的 时 间 ; 6 : 种 群 爆发 前 的 种 群 大 小 ; 9, : 种 群 爆发 后 的 种 群 大 小 ; RAG: the Hapending’s Raggedness 指 
数 ; SSD: 偏离 总 方差 ;Tajima's D, Fu's Fs; 中 性 检验 指数 ;Expansiontime (mya): 种 群 扩张 时 间 〈 百 万 年 前 ) 


T, time in number of generations elapsed since the sudden expansion episode; 6) and 0 are pre-expansion and post-expansion populations size; 








RAG, the Harpending’s Raggedness index; SSD, sum of squared deviations 


3 讨论 

本 文 的 研究 表明 青海 当归 ITS 单 倍 型 在 其 地 
理 分 布 区 存在 明显 的 谱系 地 理 格局 。Qiu 等 
(2011) 结合 前 人 的 研究 ， 总 结 出 青藏 高 原 植物 
居 群 扩张 的 模式 大 致 存在 以 下 两 种 情况 ,一 是 
植物 遵循 “冰期 生境 收缩 ， 冰 期 后 扩张 的 模式 
(Hewitt, 2000)。 青 藏 高 原 边缘 低 海 拔 的 地 区 成 
为 物种 在 冰期 的 避难 所 ， 并 保留 了 祖先 单 倍 型 ， 
而 在 青藏 高 原 台 面 上 则 只 能 检测 到 少数 几 种 单 倍 
型 ， 二 是 在 第 四 纪 冰 期 ， 植 物 在 青藏 高 原 台 面 上 
仍 有 部 分 居 群 存活 下 来 并 且 保 留 了 部 分 祖先 单 倍 
型 ， 青 藏 高 原 台 面 的 局 部 地 区 成 为 避难 所 ， 从 而 
在 高 原 面 上 有 多 种 单 倍 型 。 青 海 当 归 的 单 倍 型 分 
布 情况 更 倾向 于 第 一 种 情况 ， 青 藏 高 原 东南 部 的 
eR ( 红 原 ， 瓦 切 ， 红 星 ， 唐 克 ， 松 潘 ) 有 着 
较 高 的 遗传 变异 和 单 倍 型 特有 性 。 许 多 研究 表明 
一 个 物种 的 遗传 多 样 性 在 冰期 避难 所 与 间 冰 期 扩 
散 的 居 群 间 会 存在 差异 ， 避 难 所 的 居 群 往往 具有 
更 高 的 遗传 多 样 性 ， 而 扩散 地 区 的 居 群 遗传 多 样 
性 相对 较 低 ， 这 是 因为 当 物 种 在 间 冰 期 从 避难 所 
向 外 扩张 时 ， 伴 随 着 奠基 者 效应 ， 居 群 可 能 会 受 
到 一 系列 瓶颈 作用 的 影响 ， 从 而 使 得 迁移 后 的 
居 群 遗传 多 样 性 降低 (Avise, 2000; Petit 等 ， 
2003). 。 同 时 ，Petit 认为 在 避难 所 地 区 的 居 群 ， 
不 仅 有 高 的 遗传 变异 ， 而且 有 单 倍 型 的 特有 性 
(Petit 等 , 2003 ) 。 因 此 我 们 的 研究 结果 表明 青藏 
高 原 东南 部 的 横断 山区 很 可 能 是 青海 当归 在 第 四 
纪 冰 期 的 避难 所 。 高 原 台面 的 居 群 都 具有 古老 的 
单 倍 型 HL, 序列 变异 少 ， 一 种 可 能 是 由 于 青海 
当归 种 群 在 冰期 后 扩张 的 过 程 中 经 历 了 奠基 者 效 
应 ,现存 于 高 原 台面 的 青海 当归 种 群起 源 于 第 四 
纪 冰 期 存在 于 高 原 东南 部 避难 所 种 群 ， 种 群 共 同 






























































经 历 了 一 次 冰期 后 的 回迁 ， 在 回迁 过 程 中 有 优势 
的 单 倍 型 Hl 的 个 体 首 先 占 领 青藏 高 原 适 宜 的 生 
存 环境 ， 并 不 断 扩大 其 分 布 范围 ， 黄 基 者 效应 使 
得 所 有 的 种 群 都 拥有 相同 的 单 倍 型 HI ; 另 一 种 
可 能 是 现存 于 高 原 台面 的 这 些 种 群 是 第 四 纪 冰 期 
后 的 残余 种 群 ， 即 高 原 台 面 分 布 的 种 群 在 冰期 原 
地 保留 ， 瓶 贷 作用 造成 了 种 群 内 遗传 变异 的 减 
少 。 但 根据 遗传 漂 变 对 单 倍 型 的 随机 选择 ， 高 原 
台面 的 种 群 应 该 固定 几 种 不 同 的 单 倍 型 ， 并 且 单 
倍 型 频率 也 会 有 显著 的 不 同 ， 种 群 间 的 遗传 变异 
较 大 。 然 而 我 们 的 研究 结果 表明 ， 这 种 模式 在 高 
原 台 面 是 不 存在 的 ， 高 原 台面 几 个 种 群 间 的 遗传 
分 化 非常 小 ， 因 此 ,“ 高 原 台 面 分 布 的 青海 当归 种 
群 在 冰期 原 地 保留 ”是 不 可 能 的 。 而 青藏 高 原 东 
南部 的 地 形 为 其 充当 青海 当归 第 四 纪 冰 期 的 避难 
所 也 提供 了 有 利 条 件 ， 当 冰川 来 临时 ， 秦 岭 、 云 
贵 高 原 和 大 巴山 像 屏 障 一 样 前 弱 了 冰川 的 作用 ， 
同时 由 于 四 川 盆 地 内 陆 与 云贵 高 原 温暖 多 十 的 气 
候 的 影响 ， 西 部 南北 走向 的 河谷 以 及 大 小 盆地 为 
该 地 区 植物 迁移 和 保存 提供 了 有 利 环境 ， 也 成 为 
了 许多 植物 的 冰期 “避难 所 ”( 吴 征 锚 ，1980，; 
Yang 等 , 2012) 。 第 四 纪 冰 期 被 划分 为 4 个 阶段 : 
Xixiabangma (1.17-0.80 百 万 年 前 ) Naynayxungla 
(0.72-0.50 百 万 年 前 )、Guxiang (0.30-0.13 百 
万 年 前 ) 和 Baiyu (0.01-0.07 百 万 年 前 )， 在 第 
二 次 冰期 Naynayxungla (0.72-0. 50 百 万 年 前 ) 
以 后 ,冰川 就 变 得 不 是 很 广泛 和 强烈 了 (Zheng 
等 , 2002) ， 青 海 当 归 种 群 扩张 的 时 间 大 约 是 在 
0. 138 百 万 年 前 ， 扩 张 的 时 间 在 第 二 次 大 冰期 以 
后 ， 这 就 说 明 随 着 冰川 的 减弱 、 青 藏 高 原 气 候 的 
好 转 和 温度 的 回升 ， 青 海 当 归 扩 张 了 种 群 的 范围 。 

在 青海 当归 居 群 中 ，AMOVA 结果 显示 居 群 
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间 的 遗传 分 化 为 0. 67。 结 果 表 明 青 海 当 归 种 内 
变异 主要 来 自 于 居 群 间 的 遗传 分 化 ， 尤 其 是 青藏 
高 原 东南 部 边缘 居 群 间 的 遗传 变异 。 位 于 青藏 高 
原 东南 缘 的 松潘 和 唐 克 等 居 群 都 有 各 自 的 私有 单 
倍 型 ， 而 松潘 居 群 与 其 他 居 群 没有 任何 共享 单 倍 
型 ， 这 可 能 是 由 于 生境 的 片段 化 所 引起 的 。 青 茂 
高 原 的 边缘 存在 一 系列 的 高 山 和 峡谷 ， 它 们 直接 
影响 了 青海 当归 的 种 子 扩散 范围 ， 阻 隔 了 居 群 之 
间 的 基因 流 ， 所 以 导致 了 青海 当归 在 该 地 区 遗传 
分 化 较 大 。 松 潘 居 群 两 个 特有 单 倍 型 的 形成 有 可 
能 是 由 于 岷 山 的 隆起 成 为 了 居 群 间 基 因 交 流 的 屏 
障 ， 该 居 群 的 个 体 在 适应 高 海拔 环境 的 过 程 中 衍 
生出 了 新 的 单 倍 型 。 

综 上 所 述 ， 青 海 当 归 所 有 居 群 nrDNA ITS Fe 
列 变异 研究 的 结果 表明 ， 该 物种 冰期 的 避难 所 在 
青藏 高 原 东 南部 边缘 ,青海 当 归 种 群 的 遗传 变异 
也 主要 发 生 在 避难 所 地 区 ， 而 高 原 台 面 上 经 历 了 
共同 的 回迁 和 随后 的 黄 基 者 效应 ， 因 而 在 它 回迁 
地 区 的 居 群 拥有 共同 的 单 倍 型 HI ， 且 遗传 变 蜡 
比较 单一 。 
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